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^ , ABSTRACT 

■ 

^SJ ' Context. There are several edge-on galaxies with a known magnetic field structure in their halo. A vertical magnetic field significantly 
enhances the cosmic-ray transport from the disk into the halo. This could explain the existence of the observed radio halos. 

1 I , Aims. We observed NGC 253 that possesses one of the brightest radio halos discovered so far. Since this galaxy is not exactly edge-on 
■ (; = 78°) the disk magnetic field has to be modeled and subtracted from the observations in order to study the magnetic field in the 

r \ • halo. 

' Methods. We used radio continuum polarimetry with the VLA in D-configuration and the Effelsberg 100-m telescope. NGC 253 has 

(~ l , a very bright nuclear point-like source, so that we had to correct for instrumental polarization. We used appropriate Effelsberg beam 

^ |. patterns and developed a tailored polarization calibration to cope with the off-axis location of the nucleus in the VLA primary beams. 

I ' Observations at AA 6.2 cm and 3.6 cm were combined to calculate the RM distribution and to correct for Faraday rotation, 

y I Results. The large-scale magnetic field consists of a disk (r, (p) and halo (r, z) component. The disk component can be described as an 

, axisymmetric spiral field pointing inwards with a pitch angle of 25° ± 5° which is symmetric with respect to the plane (even parity). 

^ . This field dominates in the disk, so that the observed magnetic field orientation is disk parallel at small distances from the midplane. 

1^1 ' The halo field shows a prominent X-shape centered on the nucleus similar to that of other edge-on galaxies. We propose a model 
where the halo field lines are along a cone with an opening angle of 90° ± 30° and are pointing away from the disk in both the northern 

I , and southern halo (even parity). We can not exclude that the field points inwards in the northern halo (odd parity). The X-shaped halo 

^ ■ field follows the lobes seen in Hey and soft X-ray emission. 

1^-^ ' Conclusions. Dynamo action and a disk wind can explain the X-shaped halo field. The nuclear starburst-driven superwind may further 

QQ ' amplify and align the halo field by compression of the lobes of the expanding superbubbles. The disk wind is a promising candidate 

, for the origin of the gas in the halo and for the expulsion of small-scale helical fields as requested for efficient dynamo action. 

Key words, galaxies: individual: NGC 253 - Magnetic fields - methods: observational - galaxies: magnetic fields - galaxies: halos - 
QQ . galaxies: ISM 
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' 1. Introduction nova explosions and s tellar winds of massive stars (see e.g 



iDettmar & Soidall2006l) . 



>< 



. Direct evidence for supernova-heated galactic halos comes 

5^ . There is increasing observational evidence for the existence of f^^^ imaging and spectroscopy of diffuse soft X-ray emission. 

d ; gaseous halos around disk galaxies. They consist of different -phis has been impressively demonstrated by XMM-Newton ob- 

interstellar medium (ISM) species, mainly the diffuse ionized nervations of NGC 253 where the nuclear X-ray plume can 

gas, dust, cosmic rays (CRs), and magnetic fields (for a re- ^e explained only by star form ation withou t the contribution 

view, see e.g. | Dettmaij | 1992|) . The transport of gas from the disk ^GN (iPietsch et al.l I2OOII iBaueret al.ll2008h . Moreover, 

mto the halo has been discussed in the past in terms of galac- ^j^^jy ^^ows that the galactic starburst must drive a ther- 

tic chimneys dNorman & Ikeuchj||I989|) and superbubble blow ^.^j outflow, since there are strong indications of collision- 

outs ( iMacLow & Ferrara 1999). These models all include su- ally excited oxygen and iron L line complexes in the spectrum 

pernova explosions as the energy source that drives the forma- (|Breitschwerdtll2063l) 
tion of the halo. The effects of star formation on the ISM are 

dramatic: the heated gas in supernova remnants and accelerated In a pioneering work llpavichi d 19751) pointed out the im- 
energetic particles are injected into the base of the halo. The fun- portance of CRs for the generation of a galactic wind. In the 
damental parameter, which determines the formation of a halo, disk, the CR gas, the magnetic field, a nd the hot gas co ntribute 
is therefore the energy input by supernova explosions. Thus, the roughly the same amount of pressure (iBeck et al.lll996l) . In the 
halo does not form above the entire disk, but only at radial dis - halo, however, CRs and the magnetic field dominate. The rel- 
iances where the star formation takes place (|D^lemeraLli995|). ativistic CR gas has a larger pressure scaleheight than the hot 
Moreover, the compactness of the star-forming regions deter- gas, because its adiabatic index is 4/3 whereas it is 5/3 for the 
mines the threshold condition where the hot gas can break ou t hot g as; the magnetic field has an even larger scaleheight (Bec3 
against the pull of the gravitational field (iDahlem et al.ll20()6l) . |2007|) . Furthermore, the CR nucleons, which contain the bulk of 
This corroborates the picture of the hot gas injected by super- the energy, do not suffer from strong radiative losses like the hot 
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gas. The CR-driven galactic wind is thus another scenario which 
needs consideration when discussing gaseous halos. 

The theoretical framework for the CR transport is given 
by the combined diffusion-convection equation, which can be 
applied in the local comoving coordinate system if the bulk 
speed of the background medium is non-r elativistic ( Schlickeiser 
I2OO2I) . iBreitschwerdt et alj ( ll99lL 1 19931) applied the transport 
equation for different magnetic field configurations in the Milky 
Way. They concluded that the CR transport can be split in an en- 
tirely diffusive and an entirely convective regime, depending only 
on the local magnetic field configuration. For the lower Milky 
Way halo they assumed the magnetic field lines to be turbulently 
excited by stochastic gas motions caused by expanding and over- 
lapping supernova remnants. Thus, there will be no preferred di- 
rection of propagation of these magnetic fluctuations and hence 
no net Alfvenic drift. In the halo the magnetic field lines of the 
superbubbles begin to overlap and form "ope n" field lines that 
might be enhanced by magnetic reconnection ( lParkeiill992h . 

Various polarization studies have shown that magnetic fields 
are a very sensitive tracer of int eraction in the I SM th at is not 
visible at any other wavelength. IChvzv & Bed (12004 showed 
that the interacting pair of galaxies NGC 4038/39 (the Antennae) 
possesses a strong magnetic field of 20 juG, significantly stronger 
than non-interacting spirals. The polarized emission in many 
of the Virgo clus ter galaxies is shifted w ith respect to the op- 
tical distribution dWezgowiec et al.ll2007l) . Simulations suggest 
that this behavior can be explained by ram pressure stripping of 
galax i es moving through the intracluster medium ( ISoida et all 
I2OO6I: IVoUmer et"aLll2008l) . But a close inspection of the Virgo 
spiral NGC 4254 showed that the observed magnetic field struc- 
ture requires additional MHD mechanisms other than ram pres- 
sure stripping (Chyzv et al. 2007; Chvzy 2008). The structure of 
the magnetic field is thus an important tracer for MHD processes 
and the interaction between various ISM components, which are 
expected to be present in galaxies with winds. 

From these considerations it is clear that there is a need for 
understanding the three-dimensional structure of the magnetic 
field in the halo of spiral galaxies. That restricts the observa- 
tions to a few nearby edge-on galaxies that allow us to study 
the extra-planar magnetic field with high spatial resolution and 
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one of the brightest radio halos discovered so far, but the inclina- 
tion angle is only mildly edge-on (; - 78.5°). We assumed a dis- 
tance of 3.94 Mpc jKarachentsev et al]|2003h where an angular 
resolution of 30" corresponds to a spatial resolution of 600 pc. 
Furthermore, we used p.a. = 52° as the position angle of the 
major axis. 

This paper is the successive paper of lHeesen etall (l2009t) 
(Paper I thereafter) where the CR distribution in NGC 253 was 
found to be consistent with a vertical CR transport from the disk 
into the halo. In contradiction to this finding so far no vertical 
magnetic field has been discovered in the halo of this galaxy. The 
most deta iled discuss i on of the magnetic field structure was pre- 
sented bv lBeck et al.l (1 19941) . where a mainly disk-parallel mag- 
netic field was found in the disk and halo. This was explained by 
a strong shearing of the magnetic field due to differential rota- 
tion. The study presented in the present paper provides new sen- 
sitive observations of the large-scale magnetic field in the halo 
of NGC 253. We use the polarimetry information of the observa- 
tions presented in Paper I at different wavelengths which allow 



us to apply a correction for the Faraday rotation and thus to de- 
termine the intrinsic magnetic field orientation. 

This paper is organized as follows: we start with the descrip- 
tion of the observations and explain especially the calibration of 
the polarization in order to cope with instrumental polarization 
caused by the high dynamic range (Sect.|2]i. We present the con- 
tinuum maps of NGC 253 along with the vectors of the intrin- 
sic magnetic field and briefly describe its morphology in Sect.[3] 
Sect.m summarizes the polarization properties. In Sect.|5]we in- 
vestigate the magnetic field structure and present a model for 
the large-scale magnetic field. In Sect.|6] we discuss the conse- 
quences of our findings for the radio halo of NGC 253 and the 
observed other phases of the ISM. Finally we summarize our 
results and present the conclusions in Sect.|7] 

2. Observations and data reduction 

2.1. Effelsberg observations 

In Paper I we described our radio continuum observations of 
NGC 253 at AA 6.2 cm and 3.6 cm with the 100-m Effelsberg 
telescopeQHere we only explain the details important for the po- 
larization measurements and for the general calibration and data 
reduction we refer to Paper I. The high dynamic range (> 1000) 
due to the strong nuclear point-like source (hereafter for sim- 
plicity called the nucleus) requires several additional steps in 
the data reduction. In Paper I we explained how we removed 
the sidelobes of the nucleus via a Hogbom cleaning of the to- 
tal power maps. But the high dynamic range also influences the 
polarization via the leakage of total power emission to the po- 
larized intensity. The flux density of the so-called instrumental 
polarization is about 1.0% for the Effelsberg telescope with re- 
spect to the total power flux density tHeesea.2008[) . Thus, it can 
be neglected only for observations with a dynamic range signif- 
icantly less than 100 which is not fulfilled for our observations. 

In order to apply a correction for the instrumental polar- 
ization we used beam patterns for Stokes parameters Q and U 
obtained from deep observations of the unpolarized point-like 
source 3C84. The cleaned total power emission was convolved 
with the beam patterns and the computed instrumental contribu- 
tion was subtracted from the Stokes Q and U maps, respectively. 
Finally we obtained the map of the polarized intensity from the 
corrected Stokes Q and U maps using COMB (part of AIPS)0 
We applied a correction for the noise bias in order to preserve 
the mean zero level of the polarized intensity. The rms noise 
levels of the final maps are 200 |/Jy /beam for /16.2 cm at 144" 
resolution and lOOyuJy/beam for /13.6 cm at 84" resolution, re- 
spectivelyQ 

2.2. VLA observations 

For the polarization calibration of the VLA mosaic at A6.2 cm 
(see Paper I) we met the challenge that the nucleus is located 
at the edge of the primary beam in some of our pointings^ In 
the center of the primary beam the instrumental polarization is 



' The Effelsberg 100-m telescope is operated by the Max-Planck- 
Institut fiir Radioastronomie (MPIfR). 

- The Astronomical Image Processing System (AIPS) is distributed 
by the National Radio Astronomy Observatory (NRAO) as free soft- 
ware. 

^ All angular resolutions in this paper are referred to as the Half 
Power Beam Width (HPBW). 

The VLA (Very Large An^ay) is operated by the NRAO (National 
Radio Astronomy Observatory). 
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Fig. 1. Total power radio continuum at A6.2 cm from EfFelsberg 
observations with 144" resolution. Contours are at 3, 6, 12, 24, 
48, 96, 192, 384, 768, and 1536 x 1 mJy/beam (the rms noise 
level). The vectors show the orientation of the large-scale mag- 
netic field (not corrected for Faraday rotation). A vector length 
of 1 " is equivalent to 50 /vJy /beam polarized intensity. The circle 
in the lower left corner in this and the following figures indicates 
the size of the beam. 



less than 0.1% but at the edge of the primary beam it rises to 
more than 1 %. The variability in time due to the apparent circu- 
lar motion of the nucleus around the center in the primary beam, 
caused by the alt-azimuthal mount of the VLA antennas, com- 
plicates the situation even further. This prevents the removal of 
sidelobes with the Hogbom clean algorithm as it works only for 
non-variable sources. Thus, the standard polarization calibration 
technique applied in case of NGC 253 leads to a polarization 
map completely dominated by instrumental polarization. 

Hence, we developed a specially tailored technique for the 
pointings with the nucleus located at the edge of the primary 
beam. The polarization calibration PCAL (part of AIPS) was ap- 
plied to the nucleus itself instead of the secondary (phase) cal- 
ibrator in order to find a solution for the ofF-axis instrumental 
polarization. The appropriate beam patterns in Stokes Q and U 
allowed us to subtract the nucleus from the (u,v)-data of individ- 
ual "snapshots", which contained less than lOmin of observing 
time. This effectively removes the time-variable contribution of 
the instrumental polarization caused by the nucleus. Thus we 
were able to use the full set of (u,v)-data in order to get the best 
(u,v)-coverage which resulted in the expected sensitivity of the 
extended emission. A more det ailed descriptio n of the polariza- 
tion caUbration can be found in lHeese III (120081) . 

The pointings that are not influenced by the nucleus were 
calibrated in the standard way with the secondary calibrator us- 
ing PCAL. Inverting the (u,v)-data with natural weighting (i.e. 
Briggs Robust=8) we produced maps of all pointings in Stokes 
Q and U, which were convolved with a Gaussian to 30" reso- 
lution. The combination of the pointings was done with LTESS 
(part of AIPS) which performs a linear superposition with a cor- 
rection for the VLA primary beam attenuation using informa- 
tion from eac h pointing out to the 7 % level of the primary beam 
(lBraunll988h . The final maps are essentially noise limited with a 




00 48 15 00 47 45 3D 15 00 

RIGHT ASCENSION (J2000) 



Fig. 2. Total power radio continuum at A6.2 cm from combined 
VLA + Eff'elsberg observations with 30" resolution. Contours 
are at 3, 6, 12, 24, 48, 96, 192, 384, 768, 1536, 3077, 6144, 
12288, and 24576 x 30 //Jy/beam (the rms noise level). The vec- 
tors show the orientation of the large-scale magnetic field (not 
corrected for Faraday rotation). A vector length of I" is equiv- 
alent to 12.5 yuJy /beam polarized intensity. The spurs SI and S2 
and the point-like sources no. 2 and no. 6 are indicated. 



sensitivity of 30;uJy/beam at 30" resolution in Stokes Q and U. 
We used IMERG (part of AIPS) for the combination of the VLA 
and Effelsberg Stokes Q and U maps in order to fill in the miss- 
ing zero-spacing flux; this was done for the total power (Stokes 
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Fig. 3. Total power radio continuum at A3. 6 cm from Efifelsberg 
observations with 84" resolution. Contours are at 3, 6, 12, 24, 
48, 96, 192, 384, 768, and 1536 x 500 yL/Jy /beam (the rms noise 
level). The vectors show the orientation of the large-scale mag- 
netic field (not corrected for Faraday rotation). A vector length 
of 1" is equivalent to 50/iJy/beam polarized intensity. 
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Fig. 4. Polarized intensity at /I6.2cm from Effelsberg observa- 
tions with 144" resolution. Contours are at 3, 6, 12, 24, 48, 
and 96 times 200 /zJy /beam (the rms noise level). The vectors 
show the orientation of the Faraday corrected large-scale mag- 
netic field. A vector length of 1" is equivalent to 50/iJy/beam 
polarized intensity. 
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Fig. 6. Polarized intensity at /i3.6cm from Effelsberg observa- 
tions with 84" resolution. Contours are at 3, 6, 12, and 24 
X 100 yuJy /beam (the rms noise level). The vectors show the 
Faraday corrected orientation of the large-scale magnetic field. 
A vector length of 1" is equivalent to 50 //Jy/beam polarized in- 
tensity. 



/) map also. The map of the polarized intensity was computed 
again using POLC with a correction for the noise bias. 
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Fig. 5. Polarized intensity at /I6.2cm from combined VLA + 
Effelsberg observations with 30" resolution. Contours are at 3, 
6, 12, and 24 x 30yL(Jy/beam (the rms noise level). The vectors 
show the Faraday corrected orientation of the large-scale mag- 
netic field. A vector length of 1" is equivalent to 12.5 /vJy/beam 
polarized intensity. The extensions El and E2, the spurs SI and 
S2, the relative minima Dl and D2, and the point-like source 
no. 6 are indicated. 



3. Morphology of the polarized emission 



Table 1. Maps of NGC 253 presented in this paper of total power 
radio continuum (TP) and polarized intensity (PI). 



A 


Telescope 


Resolution 


TP rms PI rms 


[cm] 






[mJy/beam] 


6.2 


VLA -1- Efi^elsberg 


30" 


0.03 0.03 


6.2 


Effelsberg 


144" 


1.00 0.20 


3.6 


Efi^elsberg 


84" 


0.50 0.10 



In this section we discuss the morphology of the polarized 
emission and the magnetic field orientation. In the figures, the 
first contours indicate emission at three times the (apparent) rms 
noise level (2"'"' x 3 rms; m = 1, 2, 3, . . .). The polarization vec- 
tors (E-vectors) were rotated by 90° in order to display the orien- 
tation of the magnetic field vectors (B -vectors). We have clipped 
the vectors below two times the rms noise level in polarized in- 
tensity. The optical background image in Figs.[T][3]is from the 
DSS (when the radio continuum map is smaller than the back- 
ground image, its boundary is indicated by dashed lines)Q A 
summary of the polarization observations presented in this pa- 
per is given in Table [T] 

In Fig. [1] we present the distribution of the total power ra- 
dio continuum emission together with the B-vectors from the 

' The compressed files of the "Palomar Observatory - Space 
Telescope Science Institute Digital Sky Survey" of the northern sky, 
based on scans of the Second Palomar Sky Survey are ©1993-1995 
by the California Institute of Technology and are distributed herein by 
agreement. 
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Table 2. Integrated flux densities. 



A [cm] 


Instrument 


Se [mjy] 


Sp [mJy] 


F[%] 


6.2 


VLA + Eff . 


1440 ± 70 


80+12 


5.6+1.1 


6.2 


VLA 


1160 + 60 


64± 12 


5.5 ± 1.3 


6.2 


Effelsberg 


1430 ± 70 


87 ± 18 


6.1 ± 1.6 


3.6 


Effelsberg 


680 ± 40 


60 ± 10 


8.8 + 2.0 



/I6.2cm Effelsberg observations. The low resolution of 144" 
shows the magnetic field structure at the largest scales. The mag- 
netic field is disk parallel only in the disk plane. Further away 
from the disk, it has a significant vertical component. 

From the combined VLA + Effelsberg observations at 
A6.2 cm we created the map shown in Fig.|2] with 30" resolu- 
tion. At the location SI the magnetic field orientation is almost 
perpendicular to the disk. This is an example of a "radio spur" 
where the large-scale magnetic field opens from a disk parallel 
to a vertical field. At its end away from the disk the point-like 
source no. 6 i s located; this pol arized source likely is a back- 
ground AGN (ICarilli et al.lll992l) . The second radio spur S2 has 
a magnetic field orientation of 45° with respect to the major axis; 
the unpolarized point-like source no. 2 is located at its end away 
from the disk. 

We show the /I3.6cm Effelsberg map in Fig.|3] At a reso- 
lution of 84" this map reveals no new details. But it is impor- 
tant for the determination of the Faraday rotation as shown in 
Sect. 14:2] 

The distribution of the polarized intensity from the /16.2 cm 
Effelsberg observations is presented in Fig.|4] It is very different 
from the total power distribution that can be be described by a 
thick radio disk with a vertical exponential profile. Fig.|4]shows 
a thick disk formed by extensions El, E2 apparent in the VLA -H 
Effelsberg map (Fig.|5]l. 

Figure |5] shows the distribution of the polarized intensity 
from the combined A6.2 cm VLA + Effelsberg observations. The 
bulk of the polarized emission arises in the disk and extends 
into the halo. In contrast, in the outer parts of the disk the po- 
larized emission is weak at this resolution. The polarized emis- 
sion extends to large vertical heights above the inner disk. The 
Effelsberg /I3.6cm map of the polarized intensity presented in 
Fig.|6]is similar to Fig.|5]where the emission is slightly more con- 
centrated to the disk than at /16.2 cm. This can be explained by 
the smaller scaleheight of the total power emission at this shorter 
wavelength due to higher synchrotron losses that influences also 
the polarized emission (see Paper I). 

The polarized intensity has a minimum between the radio 
spur SI and the extension E2. Another minimum of polarized 
intensity is located near the radio spur S2. We will refer to these 
regions as "depolarized regions". The large-scale magnetic field 
is disk parallel near the galactic midplane whereas at some loca- 
tions we find a significant vertical component. Hence, the large- 
scale magnetic field may consist of two components, one disk 
parallel and one vertical. 

4. Polarization properties 

4.1. Degree of polarization 

Linearly polarized synchrotron emission is observed from CR 
electrons spiraling in an ordered magnetic field. No linearly po- 
larized signal is detected from isotropic turbulent magnetic fields 
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Fig. 7. Degree of polarization at A6.2 cm from combined VLA + 
Effelsberg observations with 30" resolution. Contours are at 1, 
2, 4, 8, 16, 32, and 64 x 1 %. The depolarized regions Dl and 
D2 are indicated. 

with randomly distributed directions^ Hence, we can use the po- 
larization degree as a measure of the ratio between the ordered 
and the turbulent magnetic field. 

The integrated polarized flux densities were obtained by in- 
tegration in ellipses which include the extra-planar emission as 
described in Paper I (see Table|2]i. The error of the flux densities 
was calculated as the quadratic sum of a 5 % calibration error 
and the baselevel error. The polarization degree P is the ratio of 
the flux density of the extended total power emission 5 e (nucleus 
subtracted) to that of the polarized emission Sp. The polarization 
degree is between 6 % and 9 %, which is similar to values found 
in other galaxies. 

In Fig.|7]we present the distribution of the degree of polar- 
ization at A6.2 cm. It has a minimum in the midplane where it is 
smaller than 10 %. Values larger than 20 % are found at the radio 
spur S 1 and and at the extensions El and E2. The depolarized re- 
gions Dl and D2 are prominently visible as local minima. 

4.2. Distribution of the Rotation Measure 

The polarization angle of the magnetic field vectors is changed 
by the Faraday rotation which depends on the square of the 
wavelength. If the Faraday depolarization is small, we can cal- 
culate the Rotation Measure (RM) from the rotation Ax of the 
polarization vector between two wavelengths: 

A;^^ = RM ■ (V - ^2). (1) 

In Fig. [8] we present the RM distribution between AA 6.2 cm 
and 3.6 cm at 144" resolution. We see a gradual change of posi- 
tive RMs in the northeastern half of the galaxy to negative RMs 

* Note that the observed B-vectors of linearly polarized emission 
can trace either regular magnetic fields (i.e. preserving their direction 
within the telescope beam) or anisotropic fields (i.e. with multiple field 
reversals within the beam). Anisotropic fields can be generated from 
turbulent fields by shear or winds. To distinguish between these two 
components, additional Faraday rotation data is needed. The fields ob- 
served in polarization are called "ordered" throughout this paper. 
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Fig. 8. RM distribution between AA 6.2cm and 3.6cm from 
EfFelsberg observations at 144" resolution. Contours are at 
-200, -100, 0, 100, and 200 x Iradm"^ -phe vector maps at 
both wavelengths were clipped below 4x the rms noise level in 
polarized intensity prior to the combination. The sectors for the 
azimuthal RM variation are also shown (see Sect. l5.3l l. 



in the southwestern half. Strong negative RMs (< -200 rad m ^) 
in region D2 are an a rtifact due to strong Faraday depolarization 
(ISokoloffet al.|[T998h . 

According to the RM distribution the large-scale magnetic 
field has a line-of-sight component that is pointing to the ob- 
server in the northeastern half and pointing away from the ob- 
server in the southwestern half. This agrees with a spiral mag- 
netic field in the disk as expected for a mean-field a-Q. dynamo. 
For an axisym metric mode the ma^gnetic field spirals in a uni- 
form direction (iBarvshnikova et aTlll987l) . Rotation curves show 
that the northeastern half is blueshifted whereas the southwest- 
ern half is redshifted (iPencelllQSUlPuche et al.lll991l) . Thus, the 
direct ions of the ve l ocity a nd magnetic field are opposite to each 
other iKrause et al. I (ll989h showed that for such a case the spiral 
magnetic field is pointing inwards. 

Using the RM distribution we corrected the magnetic field 
orientation for the Faraday rotation. This was also done for the 
combined VLA + Effelsberg map at A6.2 cm shown in Fig.|5] 
although we do not have a RM map with sufficient resolution. 
We corrected the polarization angles for angular scales larger 
than 84" using a combined VLA + Effelsberg RM map (not 
shown). The EfFelsberg polarization maps at AA 6.2 cm (Fig.|4]i 
and 3.6 cm (Fig.|6]l could be corrected without loss of resolution. 

The effect of the Faraday correction is especially important 
at the the radio spur SI: the magnetic field vectors are almost 
perpendicular to the disk after applying the Faraday correction (a 
RM of 100 rad m"^ corresponds to a rotation of 15°). In general 
the Faraday correction made the magnetic field orientation more 
disk parallel in locations near the disk plane (compare Figs.|2] 
and|3. 



Table 3. Parameters used for the ASS model of the disk magnetic 
field. 



Parameter 


Value 


Notes 


inclination angle ; 


78.5° 


optical disk 


pitch angle 


25° 


fitted to PI 


magnetic field FWHM 


13 kpc / 26kpc 


fitted to PI profile 


magnetic field scaleheight 


1.6 kpc 


thin radio disk 


ordered magnetic field 


3.0 yUG 


fitted to RM 


electron FWHM 


13 kpc 


fitted to Ho- profile 


electron scaleheight 


1.4 kpc 


fitted to Ha profile 


electron density 


0.05 cm-^ 


from Ho- 



5. The magnetic field structure 

5. 1 . Axisymmetric model for the disk magnetic field 

In Sect.[3]we suggested that the observed large-scale magnetic 
field is the superposition of a disk parallel and a vertical compo- 
nent. We will refer to the disk-parallel magnetic field with radial 
and azimuthal (r, 0) components as the disk magnetic field. The 
vertical magnetic field with radial and vertical (r, z) components 
we will refer to as the halo magnetic field. 

As a model for the disk magnetic field we use a spiral mag- 
netic field corresponding to an axisymmetric spiral (ASS) mode 
of a galactic dynamo, which is symmetric with respect to the 
plane (even parity). The inclination angle of the disk was pre- 
scribed as i = 78.5°, which is that of the optical disk. For a pitch 
angle of 25° + 5° the polarized intensity of our model resembles 
the ob served distribution (the optical pitch angle is 20° dPencd 
Il98ll) ). We fitted the observed polarized intensity along the ma- 
jor axis with a Gaussian distribution consisting of an inner disk 
(FWHM = 6.5 kpc) and an outer disk (FWHM = 13 kpc). For 
the vertical polarized emission profile we chose the exponen- 
tial scaleheight of the synchrotron emission from the thin disk 
of 0.4 kpc at /I6.2cm (Paper I). We assumed energy equiparti- 
tion between CRs and the magnetic field (where the CR num- 
ber density is «cr ^ B-), so that the Gaussian FWHM and the 
exponential scaleheight of the magnetic field are two and four 
times larger, respectively, than that of the synchrotron emission. 
The magnetic field has a FWHM of 13 kpc and 26 kpc in the 
inner and outer disk, respectively, and a vertical scaleheight of 
1.6 kpc. The profiles of the polarized flux density of the model 
and the observations are presented in Fig.|9l A summary of the 
model parameters can be found in Table[3] 

We produced maps of Stokes Q and U by integrating along 
the line-of-sight 

Q oc ncR^i ■ cos(2i/r) U oc hcrBi ■ sin(2i/r), (2) 

where is the component of the magnetic field perpendicular 
to the line-of-sight and tfr is the orientation angle of the local po- 
larization vector The maps were used to calculate the polarized 
emission and the orientation of the magnetic field. At this point 
we do not take any Faraday rotation into account but we will 
do so in Sect. l5.3l The ASS model for the disk magnetic field is 
shown in Fig.fTOl 

The model of the disk magnetic field clearly shows why the 
polarized emission is not symmetric with respect to the major 
and minor axis: the locations of the maxima of the polarized 
intensity are shifted due to the pitch angle in counterclockwise 
direction from the minor axis. Hence, the polarized intensity dis- 
tribution resembles an S-shape, as observed. The depolarized re- 
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Fig. 9. Profiles of the polarized intensity along the major (left) and minor axis (right). The black symbols show the measured 
intensities from combined VLA + Effelsberg observations at /I6.2cm. The solid (red) lines show the distribution derived from the 
disk magnetic field model. The dashed (blue) lines show the model of the combined disk and halo magnetic field (see Sect. 15. 31 ). 
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Fig. 10. Modeled polarized intensity at A6.2 cm of the axisym- 
metric magnetic field (ASS) of the disk at 30" resolution. 
Contours are at 3, 6, 12, and 24 x 30|/Jy/beam of polarized in- 
tensity. The vectors indicate the magnetic field orientation where 
the length of the vectors of 1" is equivalent to 12.5 /iJy /beam po- 
larized intensity. 



gions Dl and D2 in the observed map of the polarized intensity 
(Fig -Ell can be explained by the model, too. They are located 
where different components of the magnetic field occur in one 
beam and thus cancel each other. We note that the magnetic field 
orientation is mainly disk parallel as expected for a spiral mag- 
netic field. The good agreement between the model and the ob- 
served distribution in the disk justifies the choice of the model. 

The comparison of the profiles of the polarized emission be- 
tween the observations and the model in Fig. |9] shows that along 
the minor axis the halo shows up as additional emission. The 
halo magnetic field is investigated in the next section. 
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Fig. 11. Polarized intensity at A6.2 cm of the halo magnetic field 
with 30" resolution after subtraction of the disk field model. 
Contours and vectors are identical to Fig.[T0] The boxes for the 
integration of the orientation angle tfr are also shown. 



5.2. The halo magnetic field 

The good agreement in the disk region of the simple model 
shown in Fig.[T0]with the observations shows that the disk field 
dominates in the disk. We note that the two radio spurs are lo- 
cated near the two depolarized regions (Dl and D2 in Fig.|5]l, 
where the projected disk field has a minimum. Elsewhere the 
visibility of the halo field is reduced by the dominating disk field. 

The observed polarized emission is the superposition of the 
disk and halo component with Q - Qa + Qh and U - Ud + t/h- 
We subtracted the disk magnetic field model from the maps of 
Stokes Q and U of the combined VLA + Effelsberg observations 
at A6.2 cm. The polarized intensity and the orientation of the halo 
magnetic field was calculated from the maps of Qh and Uh and 
is presented in Fig.fTTI 
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The halo magnetic field clearly resembles an X-shape cen- 
tered on the nucleus. Its orientation is mirror symmetric both 
to the major and minor axis. X-shaped halo fields have been 
observed in several ed ge-on galaxies with inclination angles 
> 85°, almost edge-on jKrause et al.ll2006l:lKrausell2007l) . This 
is the first galaxy with an X-shaped field which is only mildly 
edge-on (/ - 78.5°), so that the emission from the disk and halo 
are superimposed. 

In order to quantify the halo magnetic field we determined 
the mean orientation angle i/r in the four boxes shown in Fig.fTTI 
We integrated Stokes Q and U and calculated the orientation 
angle of the magnetic field by 



200 



1 1 




- arctan 


1) 


2 1 





<.a. + 90' 



(3) 



where p.a. - 52° is the position angle of the major axis and -h90° 
was included in order to convert the electric polarization angle 
into the magnetic field orientation. We find a mean orientation 
angle of if/ - 46° + 15° with respect to the major axis. 

5.3. Modeling the Rotation Measure distribution 

An analysis of the RMs, averaged in sectors, as a function of the 
azimuthal angle gives further information about the structure 
of the large-scale magnetic field. The sector integration was ap- 
phed between a galactocentric radius of 1 .4 kpc and 6.7 kpc with 
a spacing of 20° in the azimuthal angle (see Fig. [8]). An effective 
inclination of 60° describes the observed distribution of the po- 
larized emission in the disk and halo at A6.2 cm with 144" reso- 
lution (Fig.|4]i. We chose the RM distribution at 144" resolution, 
because the low resolution favors the study of the large-scale RM 
structure. 

The crosses in Fig. [12] give the azimuthal RM variation to- 
gether with their errors. There is a broad maximum between 
270° and 60° and the minimum is near = 180°. This confirms 
the result previously obtained from the morphology of the RM 
distribution: the northeastern half of the galaxy contains positive 
RMs (-90° < < 90°) whereas the southwestern half of the 
galaxy contains negative RMs (90° <0 < 270°). 

The disk magnetic field is dominating the polarized emis- 
sion. Therefore we will use the ASS model of the disk magnetic 
field in order to make a model for the RM distribution. We inte- 
grated the polarization vector along the li ne-of-sight, so that w e 



take Faraday depth effects into account (ISokoloff et al. 19981). 
From the Emission Measure EM - lOcm'^ pc jHoopes et alj 
119961) we derived an electron density of - 0.05 cm ^, where 
we assumed a pathlength of 6.5 kpc. From the Ha distribution 
we derived a Gaussian distribution of the electron density along 
the major axis with a FWHM of 13 kpc and a vertical scaleheight 
of 1.4 kpc. The RM distribution of the ASS model is presented 
in Fig. [13] It is notably asymmetric where the minimum has a 
larger amplitude than the maximum. 

Averaging Stokes Q and U of our ASS magnetic field 
model in sectors provides the azimuthal RM variation shown 
as solid line in Fig.[T2] We found reasonable agreement with 
the observed variation for an ordered magnetic field strength 
of 3.0 + 0.5 /iG. From pulsar RMs we obtain a foreground of 
RMfor = 10 + 5radm"^ at the position of NGC 25 3, which is 
near to the Galactic south pole (iNoutsos et al.ll20d8l) . Compared 
with the amplitude of the RM variation the foreground RM is 
X! 10 %, which we further neglect in our analysis. 

The azimuthal RM variation of a spiral ASS field in an RM 
screen or a homogeneous, mildly inclined emitting layer can be 
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Fig. 12. RM variation as a function of the azimuthal angle 0. The 
crosses give the observed RMs averaged in sectors as shown in 
Fig. [8] and their errors. The solid (red) lines represents the ASS 
model of the disk magnetic field together with its error interval. 
The dashed (blue) lines shows the model of the combined disk 
and halo magnetic field (and its error interval) which are both of 
even parity. 
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Fig. 13. RM distribution of the axisymmetric spiral (ASS) of the 
disk magnetic field at 144" resolution. Contours are at -150, 
-100, -50, 0, 50, and 100 x 1 radm"^. The map was clipped at 
0.2 mJy/beam in polarized intensity. 



described by a simple cosine, with a pha se shift A0 equal t o 
the pitch angle of the magnetic field spiral (iKrause et al. I I1989I) . 
However, the strongly inclined disk of NGC 253 requires to take 
into account the distribution of Faraday depths along the hne-of- 
sight. 

Our model can explain the observations much better than the 
cosine variation of an RM screen. For the screen the amplitudes 
of the maximum and the minimum are equal, whereas our model 
can reproduce both the maximum with RM a; 1 20 rad m"^ and 
the minimum /?M x -180radm"-. Furthermore, the RM screen 
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Fig. 14. Modeled polarized intensity at A6.2 cm of the com- 
bined even disk and even halo magnetic field at 30" resolution. 
Contours are at 3, 6, 12, and 24 x 30/^Jy/beam of polarized in- 
tensity. The vectors indicate the magnetic field orientation (cor- 
rected for Faraday rotation) where the length of the vectors of 1" 
is equivalent to 1 2.5 /iJy /beam polarized intensity. This Figure 
corresponds to the model shown in Fig. IA.5l 
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Fig. 15. Model RM distribution for the combined even disk and 
even halo magnetic field at 144" resolution. Contours are at 
-150, -100, -50, 0, 50, 100 and 150 x 1 radm^l This Figure is 
identical with the Fig.lA.: 



predicts the maximum at AO but our model shows that the max- 
imum is very broad, in agreement with the observations. This 
can be explained by the thickness of the disk (not a thin emitting 
layer) and the inclination of the disk together with the smearing 
by the observation beam. The difference in amplitude between 
the maximum and minimum is a strong function of inclination. 

The model of the disk magnetic field alone fits already to the 
observed RM variation (reduced - 6.5). But so far we took 
not into account the halo magnetic field that we also observe and 
which contributes significantly to the polarized emission. With 
observations at only two wavelengths there is no direct way to 
decompose the observed RM of the disk and halo magnetic field, 
because they are overlapping along the line-of-sight. Still, we 
can produce a model including the halo and disk magnetic fields. 

For the halo magnetic field we propose a model where the 
field lines are along a cone above and below the galactic plane, 
so that in projection to the plane of the sky the field forms an X- 
shape. From the orientation angle of the halo magnetic field as 
derived in Sect. |5.2| we deduce an opening angle of the cone of 
180° - 2 X ^ = 90° + 30°. Note that the contributions to the RM 
do not cancel if the line-of-sight components of the field change 
sign along the line-of-sight, as is the case in our cone model, 
because of different pathlengths towards the cone's far and near 
sides. Hence, we see the RM mainly from the front half of the 
cone. 

For the same reason we can only see the field direction in 
the southern halo, because it lies in front of the bright disk and 
acts therefore as a Faraday screen. Our models show that the 
magnetic field points away from the disk in the southern halo. In 
case that the field points also away from the disk in the north- 
ern halo, the halo field has even parity. Otherwise, if the field 
points towards the disk in the northern halo, the halo field has 
odd parity. Because the field direction in the northern halo has 
only a small influence on the RM distribution the difference of 
the azimuthal RM variation is also small. We tested all possible 




Fig. 16. Sketch of the observable components of the large-scale 
magnetic field structure. In the disk, the ASS spiral magnetic 
field is even and pointing inwards. In the halo, the magnetic field 
is also even and pointing outwards. Dashed lines indicate com- 
ponents on the rear side. The dotted part of the halo magnetic 
field is not observed. 



combinations of disk and halo magnetic field configurations. The 
result is summarized in a catalog in AppendixlAl The disk mag- 
netic field is of even parity, where the field is pointing in the same 
direction above and below the galactic midplane. An odd disk 
magnetic field in combination with any halo magnetic field leads 
to strong local gradients in the RM distribution which are not ob- 
served. The even disk field together with an even halo field fits 
slightly better to the observations with a reduced = 4.9 (see 
Fig.[T2li than with an odd halo field with;^'^ = 13 (see Fig. |A.7| i. 
The polarized intensity and the RM distribution of our favored 
model for the combined even disk and even halo magnetic field 
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is shown in Figs. [14] and [15] A sketch of the magnetic field lines 
is presented in Fig. |16Fl 

We investigated also the possibility of an azimuthal halo 
ma gnetic field component which is predicted by some ha lo mod- 
els iDaHa^cchia^lchiyillOOilmniXAb^ It was 
included into the conical magnetic field model and the RM vari- 
ation of the combined disk and halo model was studied. It turned 
out that the azimuthal component - if there is any - must point 
in the opposite direction to that of the disk magnetic field. As an 
upper limit of the azimuthal component we estimate 60 % of the 
field strength (x^ - 3.9). We note that the shape and direction 
of the model halo magnetic field is still consistent with the ob- 
servations, i.e. is not very sensitive to the relative strength of the 
azimuthal component. 

5.4. Magnetic field strengtlis 

The magnetic field strength can be calc ulated using the energy 
equipartition assumption as described in lBeck & Krausel(l2005l) . 
As in Paper 1 we use a pathlength through the total power emis- 
sion of 6.5 kpc and a CR proton to electron ratio of K - 100, 
and a nonthermal radio spectral index of 1.0. With a typical 
nonthermal total power flux density of 8mJy/beam at /I6.2cm 
with 30"resolution we find a total magnetic field strength of 
15 + 1.5 fiG in the disk. At a polarization degree of 5 % this cor- 
responds to 4.4 + 0.5 fiG for the ordered magnetic field in the 
disk. The strength of the halo magnetic field was estimated by 
fitting profiles of the polarized emission to the observations (see 
Fig. [9]). We found 4.4 + 0.5 fiG for the halo magnetic field. 

We now compare the magnetic field strengths of the ordered 
field calculated from the equipartition and from the RM analy- 
sis. The equipartition value can be larger than that of the RM 
analysis, because the anisotropic turbulent magnetic field emits 
polarized emission but does not or only weakly contribute to 
the observed RM (iBeck et al.ll2005h . From the azimuthal RM 
variation of our modeled ASS disk field we found an ordered 
magnetic field component of 3.0 + 0.5 /iG for the disk. In the 
halo, the magnetic field strength is similar with 3.0 + 0.5 /iG. 
There may be some anisotropic field component both in the 
disk and halo because the equipartition strength of the ordered 
magnetic field is about 50 % larger than that of the RM analy- 
sis. The anisotropic field component can contribute significantly 
to the polarized e mission in galaxies as sh own for NGC 1097 
jBeck et al.ll2005l) and M 5 1 (lFletchedl2009l) . 



6. Discussion 

6.1. The superwind model 

NGC 253 is a prototypical nuclear starburst galaxy, which 
is the source of an outflow of hot X-ray emitting gas 
jStrickland et al.ll2000l: iBauer et al.ll2007l) . Spectroscopic mea- 
suremen ts of Hg-emitting gas in the southern nuclear outflow 
cone bv ISchulz & Wegned (11993) give an outflow velocity of 
390 km s ' . Apart from the nuclear outflow, huge lobes of diffuse 
X-ray emission in the halo are extending up to 8 kpc away from 
the disk. These lobes are thought to be the walls of two huge 
bubbles containing a hot gas with a low density (jPietsch et al.l 



^ The magnetic field hnes are schematically shown for the observed 
components only. We miss structures that are smaller than our resolu- 
tion element and regions where not enough CR electrons are present. 
Therefore, all field lines are closed by components that we do not ob- 
serve. 



120001: [Strickland et al.ll2002h . Sensitive XMM-Newton observa- 
tions revealed that indeed the entire bubbles are filled with X- 
ray emitting gas (Bauer et al. 2008). Similar structures are found 
in other galaxies with intense star formation, with M 82 as the 
most prominent exam ple, and are now known as superwinds 
jHeckman et al.ll2000l) . 

Figure[T7]shows the difiiise Ha emission in greyscale. There 
are several Ho- plumes extending from the disk into the halo. We 
note that the large plume east of the nucleus corresponds to the 
radio spur SI (see Fig.[5]l. The diffuse X-ray emission in Fig. [18] 
shows a similar structure with the lobe at the extension El ex- 
tending furthest into the halo. The abundance of heated gas has 
an asymmetry with respect to the minor axis: the northeastern 
half contains significantly more diffiise Ha and soft X-ray emis- 
sion than the southwestern half which cannot be explained by a 
symmetri c superwind. This asy mmetry is also observed in H I 
emission (iBoomsma et al.ll2005l) which surrounds the superbub- 
bles. An explanation is attempted in Sect. 16.41 

6.2. The disk wind model 



The halo magnetic field allows CRs to stream along the field 
lines from the disk into the halo. In Paper 1 we determined the 
vertical CR bulk speed as vcr = 300 + 30kms"' which is sur- 
prisingly constant over the entire extent of the disk. This shows 
the existence of a disk wind in NGC 253. CRs cannot stream 
faster than with Alfven speed with respect to the magnetic field, 
because they are scattered at s elf-excited Alfven waves via the 
so-called streaming instability dKulsrud & Pearcdl 19691) . A typ- 
ical magnetic field strength of B » 15 ^G and a density of the 
warm gas of « ^ 0.05 cm"^ (Sect. l5.3l l leads to an Alfven speed 
of L-A = Bl K 150 km s"'. The super- Alfvenic CR bulk 

speed requires that the CRs and the magnetic field are advec- 
tively transported together in the disk wind. CRs stream with 
Alfven speed with respect to the magnetic field that is frozen into 
the thermal gas of the wind. The measured CR bulk speed is the 
superposition vqr - ww + of the outflow speed of the thermal 
gas f w in the wind and the Alfven speed i/a (iBreitschwerdt et all 
I2OO2I) . Adopting this picture the wind speed at the galactic mid- 
plane in NGC 253 is u„ ^ 150 km s '. 

The Alfven speed is difficult to estimate, so that the uncer- 
tainty in the actual value of the wind speed is large as well. 
Furthermore, we could measure only the average wind speed 
between 1 kpc and 5 kpc above the plane. Therefore, the bulk 
speed near the galactic midplane may be lo wer than 150 km s~ ' . 
Studies of CR transport in a galactic wind bv lBreitschwerdt et alj 
(119911) and follow-up works assumed a wind speed of only 
a few lOkms ' at the starting point in the galactic midplane 
which acce lerates to the escape velocity at larger distances 
to the disk dBreitschwerdt et alj|l993l: IZirakashvili et al.l 1 19961: 
IBreitschwerdt et al.l |2002|) . On the oflier hand. [Everett et'all 
(I2OO8I) found a high wind speed of 170 km s ' in the plane of 
the Milky Way in order to explain the soft X-ray emission in the 
halo. Our value measured in NGC 253 does not exclude any of 
these models. It will be left to future observations to provide an 
accurate estimate of the wind speed at different heights above 
the galactic plane. 

If the large-scale magnetic field is frozen into the thermal gas 
in the wind the halo magnetic field has an azimuthal component 
and winds up in a spiral. The field configuration is similar to that 
of the sun with the P arker spiral magnetic field (fParker 1958|; 
IWeber & Davis|[T967h . The azimuthal component increases with 
height, so that above a certain height the magnetic field is almost 
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purely azimuthal. Up to the Alfvenic point the stiff magnetic 
field lines corotate with the underlying disk. If the azimuthal 
component of the halo magnetic field is small, the observable 
radio halo is within the corotating regime and the Alfven radius 
is larger than the vertical extent of the halo. As an azimuthal 
halo component cannot be excluded (Sect. l5.3l l. the Alfvenic ra- 
dius may be smaller than 2 kpc, which is the scaleheight of the 
polarized emission in the halo. 

The existence of a disk wind has important consequences 
for the transport of angular momentum. As we have shown in 
Paper I the CR bulk speed is fairly constant over the extent of 
the disk, i.e. it does not depend on the galactocentric radius. 
If this is also true for the speed of the disk wind, angular mo- 
mentum is effectively tr ansported to large galactocentric radii 
dZirakashviU et al.lll996h . Moreover, the angular momentum of 
the wind per uni t wind mass is proportio nal to the Alfven ra- 
dius as shown bv lZirakashvili et akl {[996). An Alfven radius of 
«; 2 kpc agrees with their models for which they found a signifi- 
cant loss of angular momentum over the lifetime of a galaxy. We 
note that the disk wind can account for larger angular momentum 
losses than the superwind in the center (at small galactocentric 
radii). 
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Fig. 17. Halo magnetic field overlaid onto diffuse Ha emission 
as greyscale. Contours are at 3, 6, 12, and 24 x 30^(Jy/beam 
(the rms noise level). A vector length of 1" is equiva- 
lent to 12.5 /uJy /be am polarized intensity. The Ho- map from 
iHoopes et al.l (Il996h was convolved with a Gaussian to 20" res- 
olution in order to show the weak diffuse emission in the halo. 



6.3. The origin of X-siiaped ha/o magnetic fieids 

The distribution of the halo magnetic field is X-shaped in both 
orientation and intensity. As the intensity of the polarized emis- 
sion depends on the perpendicular component of the ordered 
field, a possible explanation of the intensity distribution is limb 
brightening in the conical halo magnetic field. We modeled 
this effect using a halo magnetic field lying on a cone with an 
opening angle of 90° and the inclination angle of NGC 253. 
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Fig. 18. Halo magnetic field overlaid onto diffuse X-ray emis- 
sion. The contour is at 3 x 30 //Jy /beam (the rms noise level). 
A vector length of 1" is equivalent to 12.5 //Jy/beam polarized 
intensity. The X-ray map is fro m XMM-Newton observations in 
the energy band 0.5 - 1 .0 keV jBauer et al.ll2008l) . 




Fig. 19. Halo structu re of NGC 253. Reproduced from 
iBoomsma et al.l ( |2005|) and extended. The superbubble, filled 
with soft X-rays emitting gas, expands into the surrounding 
medium (indicated by dotted lines with arrows). The halo 
magnetic field is aligned with the walls of the superbubble. 
Dashed lines denote components that are not (or only weakly) 
detected in the southwestern half of NGC 253. 

The model can reproduce indeed an increased vertical extension 
along the limbs. 

The halo magnetic field follows the lobes of the heated gas 
that are regarded as walls of huge bubbles expanding into the sur- 
rounding medium. The large-scale magnetic field may be com- 
pressed and aligned by shock waves in the walls. These shock 
waves are also able to heat a pre-existing cold halo gas, so that 
the gas becomes visible as Ha and soft X-ray emission. The cold 
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halo gas seen in HI emission surrounds the superbubbles and 
shows the same asymmetric distribution as expected if this gas 
is the source for the heated gas in the halo. A cartoon of the 
halo structure including the large-scale magnetic field is shown 
in Fig.[l9] 

Both the disk wind model and the superwind model, in con- 
junction with a large-scale dynamo action, can explain the X- 
shaped halo magnetic field. We should note, however, that such 
magnetic field structures are observed in several edge-on galax- 
ies. These gala xies show a very different level of star formation 
(lKrausell2008h . Most of them do not possess a starburst in the 
center and hence they have no superwind. While the classical a- 
O dynamo can not explain such an X-shaped halo field, model 
calculations of a galactic disk within a spherical halo including 
a galactic wind showed a similar magnetic field configuration 
dBrandenburg et all Il993l) . Hydrodynamical simulations show 
that the wind in spiral galaxies has a significant radial compo- 
nent due to the radial gradient of the gravitational potential. The 
wind flow reveals an X-shape w hen the galaxy is seen edge-on 
dPalla Vecchia & Schav3l2008l) . New MHD simulations of disk 
galaxies including a galactic wind are in progress and may ex- 
plain the halo X-shaped field dGressel et al.ll2008t iHanasz et al.l 
l2009aiibi) . The first global galactic-scale MHD simulations of a 
CR-driven dynamo give very promising results showing directly 
fliat magnetic flux is transported from the disk into the halo 
dHan asz et al. 2009cl)- Therefore we suggest that in NGC 253 the 
X-shaped halo field is connected rather to the disk wind than to 
the superwind. The superwind gas flow may be collimated by 
the halo magnetic field. 

The transport of CRs and magnetic fields has important con- 
sequences for the possibility of a working dynamo in NGC 253. 
The a-Q. dynamo relies on the differential rotation of the galactic 
disk (Q-effect) and the cyclonic motions of the ionized gas {a- 
effect). The latter one can be generated by subsequent supernova 
explosions dFerrierdl 19921) . The amplification of magnetic fields 
by an a-Q dyn amo requires expulsion of small-sc a le helical 
fields (see e.g. [Brandenburg & Subramanianll2005b . ISur et al.l 
d2007h showed that galactic winds can transport the helicity flux 
by advection. A galactic wind may be thus vital for a working 
dynamo. 

6.4. The origin of the gas in the halo 

The cold gas in the halo of NGC 253 is necessary to explain the 
halo structure with the superw i nd mo del. Its origin, however, is 
yet unknown. iBoomsma et al.l d2005h discussed a gas infall in 
NGC 253 by a minor merger to explain the asymmetric H 1 halo 
distribution with respect to the minor axis. But they noted that 
a minor merger event is unlikely because the distribution of H 1 
is very symmetric with respect to the major axis. Could the disk 
wind be the origin of the cold gas in the halo? The disk wind is 
asymmetric with respect to the minor axis. The convective north- 
eastern half indicates a strong disk wind which transports gas 
from the disk into the halo while in the diffusive southwestern 
half the wind is weaker (Paper I). 

We can only speculate about the reason why the two halo 
parts are so different. Star formation in the disk most likely plays 
an important role. The northeastern spiral arm contains signifi- 
cantly more H 1 1 regions than the sou thwestern one as is visible 
from Fig. 4b in iHoopes et al.l (Il996l) . Moreover, the amount of 
total radio continuum emission indicates a higher star-formation 
rate in this part of the disk. Since the energy input of star forma- 
tion is essential for the injection of thermal and CR gas, the disk 
wind can be more easily driven in the northeastern half. Thus, 



the disk wind is a good candidate for the origin of the cold gas 
in the halo. 



7. Summary and conclusions 

The three-dimensional magnetic field structure can be investi- 
gated using sensitive radio continuum polarimetry. Our main re- 
sults are; 

1. A disk-parallel magnetic field exists along the midplane of 
the disk. The magnetic field lines are slowly opening further 
away from the midplane. The vertical magnetic field compo- 
nent is most prominent at the edge of the inner disk where 
we find two "radio spurs", one previously known east of the 
nucleus and one newly discovered west of the nucleus. The 
magnetic field configuration can be described as an X-shaped 
pattern, as in other edge-on galaxies. 

2. The large-scale magnetic field can be decomposed into a disk 
(r, (p) and halo (r, z) component. 

3. The disk magnetic field of NGC 253 can be described by 
an axisymmetric spiral (ASS) magnetic field with a constant 
pitch angle of 25° +5° which is symmetric with respect to the 
plane. This model shows a high resemblance to the observed 
magnetic field in the disk. 

4. The distribution of the polarized intensity and the orientation 
of the halo magnetic field shows a distinct X-shaped pattern 
centered on the nucleus. Our model of the disk magnetic field 
was subtracted to construct a map of the halo magnetic field. 
We propose a model where the halo field lines are along a 
cone with an opening angle of 90° + 30° and are pointing 
away from the disk (even parity). An odd halo magnetic field 
is also possible, because we can not reliably determine the 
magnetic field direction in the northern halo. 

5. The distribution of the halo magnetic field coincides in shape 
with the extra-planar heated gas traced by Ha and soft X- 
ray emission. Possible explanations are limb brightening and 
compression of the halo magnetic field in the walls of the 
expanding superbubbles. 

6. A disk wind plus dynamo action is a promising scenario for 
the origin of the gas in the halo and for the expulsion of 
small-scale helical fields as requested for efficient dynamo 
action. The disk wind can also account for large angular mo- 
mentum and magnetic flux losses over galactic time scales. 

Future investigations should include polarimetry at low frequen- 
cies, to take advantage of studying the magnetic structure far 
away from the disk illuminated by an aged population of CR 
electrons. This should become feasible with the upcoming class 
of telescopes operating at low frequencie s like LOFAR and the 
low-f r equency SKA arra y. RM Synthesis dBrentiens & de Bruvnl 
I2005t iHeald et al.l I2OO9I) applied to multi-channel polarization 
observations is required to separate the RM contributions from 
the disk and the halo of strongly inclined galaxies. 
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Appendix A: Magnetic field models 

In order to find the best model for the magnetic field in NGC 253 
we studied the expected polarized intensities and RMs of vari- 
ous possibilities of the disk and halo magnetic field alone and 
combinations of them. The setup of the models is described in 
Sect. lS.ll with the parameters of Table |3] In Figs. lA.lllA.12l we 
present the polarized intensity and the magnetic field orientation 
(not corrected for Faraday rotation) for AA 6.2 cm (a) and 3.6 cm 
(b) at 30" resolution. Note that the Faraday corrected magnetic 
field orientation is identical for all models consisting both of the 
disk and halo field (Figs. |A.5llA.12t . The difference in the mod- 
els is only in the direction of the magnetic field that becomes 
visible in the RM distribution (c) shown with 144" resolution. 
Sketches (d) show the direction of the magnetic field in the disk 
and in the halo. Here, a "■" indicates that the field points to the 
observer and a "-1-" denotes a field pointing away. 

Figures [ATI and [A2\ show the disk magnetic field alone for 
the even and odd case. They are similar but the odd field has 
a smaller amplitude in the RM distribution. This can be under- 
stood as the scaleheight of the disk is smaller than the projected 
minor axis. We see therefore only one side of the disk (the south- 
ern one). Figs. lA.3 l and IA.4I show the halo magnetic field alone 
for the even and odd cases. Note the asymmetry of both the po- 
larized intensity and the RM. Because the magnetic field orienta- 
tion changes along the line-of-sight, there is some depolarization 
as the magnetic field vectors are not parallel to each other If the 
Faraday rotation has the same sense as the rotation of the mag- 
netic field, the depolarization is stronger. If they have opposite 
senses, the depolarization is weaker This can explain the asym- 
metry although the magnetic field is axisymmetric. 

Figures IA.5HA.8I show the models with the even disk mag- 
netic field. Among them Figs. |A.5| and |A.7| are the models whose 
RM distribution agrees best with the observations. The two mod- 
els only differ in the direction of the northern halo field. The even 
halo magnetic field (Fig. IA.5l l is our best-fit model. Figs. |A. 91 - 
IA.12I show the models for an odd disk magnetic field. They all 
have strong gradients in the RM which are not observed. 
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Fig. A.l. Even disk magnetic field. For Figs. lA.lllA.12l a) polarized flux density and magnetic field orientation at A6.2 cm with 30" 
resolution, b) likewise but for /I3.6cm. c) RM distribution between AA 6.2 cm and 3.6 cm with 144" resolution, d) sketch of the 
magnetic field direction. 




Fig. A.2. Odd disk magnetic field. 




Fig. A.3. Even halo magnetic field. The halo field points away from the disk. 
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Fig. A.4. Odd halo magnetic field. The halo field points away from the disk in the southern halo. 




Fig. A.5. Even disk magnetic field and even halo magnetic field. The halo field points away from the disk. 




Fig. A.6. Even disk magnetic field and even halo magnetic field. The halo field points towards the disk. 




Fig. A.7. Even disk magnetic field and odd halo magnetic field. The halo field points away from the disk in the southern halo. 
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Fig. A.8. Even disk magnetic field and odd halo magnetic field. The halo field points towards the disk in the southern halo. 
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Fig. A.9. Odd disk magnetic field and even halo magnetic field. The halo field points away from the disk. 
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Fig. A.IO. Odd disk magnetic field and even halo magnetic field. The halo field points towards the disk. 
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Fig. A.ll. Odd disk magnetic field and odd halo magnetic field. The halo field points away from the disk in the southern halo. 
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Fig. A.12. Odd disk magnetic field and odd halo magnetic field. The halo field points towards the disk in the southern halo. 



